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The goal of this paper is to compute the zeta function determinant for the massive Laplacian on 
Riemann caps (or spherical suspensions) . These manifolds are defined as compact and boundaryless 
_D— dimensional manifolds deformed by a singular Riemannian structure. The deformed spheres, 
considered previously in the literature, belong to this class. After presenting the geometry and 
discussing the spectrum of the Laplacian, we illustrate a method to compute its zeta regularized 
determinant. The special case of the deformed sphere is recovered as a limit of our general formulas. 



I. INTRODUCTION 



Zeta functions are a topic of ubiquitous interest due to their wide range of applicability. The basic reason for their 
popularity is that systems with boundaries or finite size effects (in general boundary conditions) can be conveniently 
studied by constructing and analyzing suitable functions of the spectrum of some operator. For instance, many 
fundamental properties of quantum fields are encoded in the effective action, that contains information about the 
(-H ] induced quantum energy-momentum tensor that is the fundamental object necessary to discuss, for instance, how the 
classical equations of motion are modified by quantum effects, or phenomena like particle creation in curved space. 
The effective action F can be expressed as 



F-lndetO, (1) 

with O being (typically) a positive elliptic differential operator. When, the operator O acts over a compact manifold, 
• the spectrum is also discrete. The effective action is then formally written as sum over the eigenvalues (the sum is, 
\^ ' in fact, an integral with measure depending on the spectral density of the eigenvalues in the continuous case) 

o 

q; F^^lnA. (2) 

o 

Written as above, the effective action is divergent and it does not have a precise meaning unless it is regularized 
^-H by using some appropriate analytic continuation procedure. One such procedure, introduced in [l-iSj, uses the zeta 
IL" ' function 

to define the effective action as the analytic continuation of the derivative C'{s) at s = 0. This is, in essense, the idea 
behind zeta function regularization and the standard textbook expression for the effective action is [4|] : 

^ = -^C'(0)-ic(0)lnM^ (4) 

where is a renormalization scale. Sometimes the above functional determinant is also addressed as zeta regularized 
determinant. 

Finding the zeta function and performing the analytic continuation is, most of the times, a non-trivial task, and 
only a few general methods are available. When the spectrum can be found explicitly (this is typically possible when 
the background manifold is highly symmetric and the differential operator in question is conformally covariant), aside 
from trivial cases, different ways have been adopted to compute the zeta function and functional determinant. A 
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classical example is that of the Laplacian on a spherical domain analyzed in ref. [H. In more general set-ups, when 
eigenvalues are known only implicitly, one may use the method developed in Q- 

Up to now, explicit results for the zeta determinant of the Laplacian have been obtained in a variety of cases with 
different methods that, aside of those mentioned of spheres (Hi and orbifolded spheres @, include projective spaces 
[SJ , balls [7|] , hyperbolic manifolds [9] , cones 0, [13] , intervals [ll| . 

One interesting case, beyond those that have been analyzed so far, is that of the deformed spheres, the simplest 
example of such geometries being a spherical cap. The deformed spheres can also be seen as generalization of lunes and 
orbifolded spheres. These type of geometries occur frequently in cosmology, since sectors of de Sitter space belong to 
this class. In general, a deformed sphere is defined as a Z?— dimensional sphere equipped with a singular Riemannian 
structure. The metric is written as 

dsl = de^ + sin^ edsl,_i , (5) 

where ds-jy_i is the metric of a (D — 1)— dimensional sphere. Ref. [T^] describes the geometry of deformed sphere and 
we address the reader to that reference for details. The zeta determinant of the Laplacian on the above manifolds has 
been recently analyzed in Ref. that introduced a general method to deal with some classes of simple and double 
abstract zeta functions allowing to compute the zeta invariants C(0) and C'(0): ^^'^ presented explicit formulas in 2 
and 3 dimensions. Alternative computations include those of ref. [5| that considered the case of spherical domains, 
spherical caps and lunes and computed the functional determinant in a variety of dimensions. Some related work in 
physics, aiming at computing the effective action in the context of quantum cosmology discussed the analogous 
problem of computing the zeta determinant for a sector of four dimensional de Sitter space. In brane world cosmology, 
ref. dealt with the computation of the zeta determinant in 5— dimensional anti-de Sitter space bounded by de 
sitter branes and in sectors of higher dimensional de Sitter [T5j . Finally, the Casimir energy for two-dimensional 
deformed spheres has been considered in ref. [iGj . 

In this work, we will focus our attention on some generalizations of the deformed spheres in which the base manifold 
^ is smooth, compact and without boundary, but otherwise unspecified. The base manifold is deformed by a singular 
Riemannian structure analogous to that of the deformed spheres and the local geometrical properties of the spherical 
suspension can be described by use of the following metric 

ds^ ^ de^ + sin^ edT.^ , (6) 

where dY,^ represents the line element of the base manifold The angular coordinate 6 is restricted to the interval 
[O,0o] with Oq g (0,7r). We will indicate this manifold as and use the terminology of Riemann Cap. Another 
commonly used name is spherical suspension. 

The goal of the paper is to present a method to compute the functional determinant for the Laplacian on f2. 
By confronting the problem with that studied in ref. [Toj . one immediately notices that one needs to consider the 
asymptotic behavior of new eigensolutions. This is not a merely technical complication and, as explicit computations 
show, the method developed in pi| does not easily generalize to this case, although, in principle, directly applicable. 
Here, we take a different route, analogous to that developed in ref. [13], that makes use of the direct analytic 
continuation at s = of the integral representation of the zeta function. The method is valid for only when the 
(uniform) asymptotic expansion of the eigenfunction satisfies a certain form, but it is general enough to encompass 
the cases treated by the method of [10] and also those of the spherical suspension. One of the aims of the paper is 
precisely to refine the results of ref. [13] in order to be able to treat the case of Riemann caps. 

This paper is divided in four main sections. Sec. |TT]is devoted to describe spectrum, eigenvalues and eigenfunctions, 
of the Laplacian on ft. In sec. IIIII we present the general method used to compute the functional determinant. The 
general results for the regularized zeta determinant of the Laplacian on are reported in sec. IIVI The case of the 
deformed spheres, considered in jl^ . is discussed in sec. IVland previous results recovered from our general formulas. 
The final section is devoted to summarize and discuss the results of our work. 

II. GEOMETRY AND SPECTRUM 



The Laplacian on f2 can be explicitly written, in spherical coordinates, as 

= — + dcot 0— + — 2 „ 
06'^ dO sm^ 9 

with Vs being the Laplacian on the base manifold The eigenvalue equation is 

52 . „ a 1 



dcot9— + -r^V^ , (7) 



dcot 6— -I Ve 4- TO 



sm 



<^„,k(^^,^») =<k'Pn,k(e,X^) , (8) 
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where the coordinates on jV are labehed by X^, with i = 1, • • • The indices n and k have been introduced 
to distinguish between the angular direction and the base manifold and are, for the moment, unconstrained. The 
quantity m is constant and represents a mass term once equation ([8]) is interpreted as the Klein-Gordon equation for 
a massive scalar field propagating on f2. We assume here that a complete and orthonormal set of of solutions to the 
eigenvalue equation on .vV exists, but leave the specific properties of the base manifold and the boundary conditions 
unspecified and carry out our analysis in general. When the base manifold is a sphere, we reproduce, using a different 
method, the results of refs. [1, [H,[l5|. A general solution to (|5]) can be written as 

<^„,k(e,X») = 0„,kWHk(X,) , (9) 
where HkiXi) are the harmonics on with degeneracy (i(k) and eigenvalues Ak satisfying the following equation 



(10) 



with > 0. Some simple manipulations allow us to recast eq. 



dcoH 



1 



de sin^ e 



J as follows 



By means of the following transformation. 



0„,kW = sin«0V„,k(0) 



(11) 



with q = {\ — d)/2, the previous equation can be cast in standard Legendre form 



d 

■^^nM{d)+COie — ^nM{^) 



2 2 

+ m — a„ 



{d-l) 



^nMiO) = . 



Introducing the quantities 



fJ-k 

VnM 



+ A2 , 



k , 



(12) 
(13) 



with (j^ = 



and noticing that [/^k] > and 5R [i^„,k] > ~l/2, two linearly independent solutions to the 



above equation can be written in terms of (see Theorem 12.1, ref. il7| ) 



(1 - .^)^/^p:|%^„., iz) 



(1 



(^) 



(14) 



we have used the symmetry property Pj^^ ~ P^i-iy defined z — cos 9. The two functions above correspond to 
the Ferrers representation of the Legendre function (see Chapt. 5, Sec. 15 of ref. [l3|) and satisfactorily cover the 
interval —l<z< 1, with cuts from (— oo, —1] and [1, -|-cxd). By looking at the z — ^ 1 behavior of the above functions, 
it is easy to see that the solution proportional to Q diverges in this limit. Thus, requirement of regularity at — 
restricts the general solution to be 



^„,k(^,x,) - (1 - zy/'p-_'ii,^^, (z)Hk(x,) 



(15) 



The solution (|15p . aside for the requirement of regularity at = 0, is general. In the following discussion, for 
definiteness, we will consider the case of Dirichlet boundary conditions, thus imposing on the eigenfunctions the 
relation 



Fk (ar. 



(16) 



Solutions to the previous equation implicitly determine the coefficients iuj„^k and, via (|13p. the eigenvalues a„^k- It 
should be easy to see that other boundary conditions can be easily be treated along the same lines as for Dirichlet. 
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III. ZETA FUNCTION ON Q 



In this section we wish to present a formahsm to compute the ^eta function for the Laplacian on f2. The method we 
adopt is different from the one described, for example, in ref. [lOj. The common step is the use of a contour integral 
representation for the spectral sum, 

Cis) = J2 di^)^n'u = E dik){uJnM + O^)-^ , (17) 
n,k n,k 

but it will differ in the method of ref. otherwise. The goal is to express the above zeta function in terms of the 
auxiliary zeta function defined on the base manifold ,A' in the spirit of ref. ^] , 

C^{s) = Y,d{^.)^^^^^ . (18) 

k 

The use of the residue theorem allows us to express the zeta function (|17p as a contour integral: 

= i E / dzi^' + '^'r'§; InFk(^) , (19) 

k 

where the circuit 7 (see fig.[T]) encloses all the roots of ([16]). As it stands, the above integral representation is valid in 
the region 3? [s] > (d + l)/2. By appropriately deforming the contour of integration along the imaginary axis and by 
utilizing the property = P^i_^, the above expression for the zeta function can be recast in the following form 



1/2+w 



(cos «o) 



C(.) = ^E'^(k) fdwiw^-ar^-^lr^PZT; 

+ J2dik)[ ±,iz^ + a')-s^lnF^(z), (20) 

where the contour is deformed as illustrated in Fig. [1] and Cr is a small rectangle enclosing the cut on the imaginary 
axis. The circuit is closed at infinity by a semi-circle, not shown in the figure, whose contribution to the contour 
integral vanishes due to the behavior of the integrand at large distance (the reader can check this behavior in the 
asymptotic expansion shown below). A simple calculation shows that the integral over Cr does not contribute to Ci^)- 
The first integral in the above expression comes, instead, from the portions C+ and C_ of the deformed circuit, and, 
by performing the change of variables f^^u^ — vo^ — cr"^ in the expression (|20|) , one readily obtains 

^ s™^ rdud ^^(eos^o) . (21) 

We want to show, now, that the previous expression is convergent when the parameter s is large enough. To this 
aim, let us use the large /i^ expansion of the generalized Legendre function [ij, [l^ [23| . Using a WKB form for the 
solution of the Legendre differential equation, one obtains a uniform asymptotic expansion for the Legendre functions 
of large order. Explicitly, one gets 



, (cos^o) ~ xl^—e''^^fiuu)-''^Y^^ , (22) 



where we have defined 



t = t{u) = — -^^^=^^^= , 1/ = tcost 
Vl + sin^ 6*0 



t{u) = 1 + In 



u sm Hq 



\/l + u'^ sin^ 6*0 + cos / 
The polynomials v4„(z^) can be found to satisfy the following relation [l| 



, 1 \ ^ 1 

arctan — — arctan f 



u J \ut cos 0Q 



(23) 



00 I 00 j-j I 00 

E A^k"A.(-) = exp - 5] A (-) , (24) 

«=o j=i 2j(2j + l)/^k ,=0 
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Figure 1: The figure illustrates the contour of integration 7 and the deformation used in the calculation. 



where B2j are the BernouUi numbers and the polynomials An{v) obey the recurrence relation 

„,2 



2(u^ + 1) 



8(u2 



1 



dv' ( ^V'^ + — - 1 - 2(^ \ 2 L -^"(^ ) 



with yio('^) = 1- For the interested reader, few of the first Aniy) are listed in the appendix of [19| and higher order 
ones can be obtained with a simple computer program. It is, a this point, straightforward to obtain an asymptotic 
expansion for the logarithm of the Legendre function. From the previous results we have 



InP" 



where the polynomials a„(:^) are defined through the cumulant expansion 



an{v) 



(25) 



y — 

^ 2j(2j - 



2j 



- + log< i + ^Mk"A, 



n=l 



E 

fc=i 



ak{v) 



(26) 



The previous expansion is uniform and valid for all ranges of variation of u with the coefficients being bounded 
functions in u. Substituting the previous expression in formula ([21]) one obtains the following asymptotic formula for 
the zeta function: 



as) 



sinTTs ■^-^ d{Vi) 



du d 1 
I 2 



In- 



t 



27r^i, 



ln(/^kw)^'' + MkT + ^ 



a„iu) 



1 /^k 



(27) 



Since the /ik are the eigenvalues of a second order partial differential operator on a smooth compact manifold, ^/K, 
they behave, asymptotically, according to the Weyl's estimate, as [2l| 



,/2 2^-V'^2dr(d/2)„ , 
Mk — lkl • 



Vol 



(28) 



This means that the spectral sum ([TS]) converges for [s] > d/2. Therefore, by using the fact that the coefficients 
a„(u) are bounded functions, proving that the above expression is convergent, for 3? [s] large enough, requires only 
trivial steps. In the convergency region for (|27p . we can safely invert the summation and integration operation and 
write: 



, . sinTTs f°° du dQ , . 



(29) 
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where for notational convenience we have defined 

giu,s) = E^(k)M.-"lnP:- ^^.^^(cos^o) • (30) 

k 

It should be clear, by looking at formula (l?71) . that in the asymptotic region defined by the expansion (|25l) . the function 
s) can be expressed in terms of the base zeta function, Cy^'{■^)■ Further properties of g{u, s) can be obtained by 
using the Mittag-Leffler expansion for the base zeta function 

1 °° C 

where f{s) is an entire function and C'p/2 are the integrated heat-kernel coefficients of the Laplacian Vs on the base 
manifold Using the above expression, one cas see that Q will develop, at most, a pole in s in the process of 
analytic continuation to s 0, i.e.: 

g{u, s) = -Gpiu) + GRiu) + 0(s) , (32) 
s 

where we indicated with Qp{u) and Gr{u) the polar and regular part of Q{u, s) for small s. 

At this point we need to perform the analytic continuation of the zeta function to s = 0. For this purpose we will 
use the following result Ts^ : 

Lemma. Let f{x) be a function defined for x > e with e > and analytic at x — e. Assume that f{x) has the 
following general asymptotic behavior for x — )■ oo; 

Ph<N 

fi^) = E (/^ + /'^l^^) + [/liogln^ + [/]-g + O(x-i) , pfc > , (33) 
fc=i 

where the .subscripts log and reg refer to the solely logarithmic and regular (non-singular) parts of f{x) in the large x 
limit. Then, there exists the analytic continuation of the integral 



i; 



where [/]- = - /(e). 



Taking the asymptotic expansion ()25p and Taylor expanding all the terms for large u, one can immediately see that 
the behavior ([55)1 is exactly reproduced. In other words, both functions Gp and Gp admit an expansion of the type 
((33)) for large u. Applying the above lemma to Gp and Gr, we obtain an expression for ^(s) and C'{s) in terms of Gp 
and Gr. More precisely for s — ^ one has 

C(^*) = -^[Gp]log + [GR]log+{[Gp]reg^Gp{0)) 

+ s [{[GRleg- Gr{0)) - duHnu^^Gp{u')^ +0{s^) , (35) 

and also 

C'(s) - [Gp],,, + I {[GrI,, GpiO)) £ du' inu'^Gpiu')^ + 0(s) . (36) 

It is well known [13, H^-H^i that for singular problems the heat kernel asymptotic expansion presents, in general, 
logarithmic terms. This non-standard behavior yields to the appearance of a pole in the spectral zeta function for 
s = 0. In our case, this singular term is proportional to [Gp]iog- As we will see in the next section, ^(s) on Ai has a 
pole at s = provided that C/k(s) has a pole at s = ~l/2. By recalling the formula (I5T|) it is straightforward to see 
that the residue of C^{s) at s = —1/2 is proportional to the heat kernel coefficient C(d+i)/2- This coefficient vanishes 
identically only if the base manifold ^ is even-dimensional and without boundary. In all other cases it is, in general, 
non zero. In these singular situations the zeta function cannot be directly used, in the standard fashion, in order to 
compute the C- regularized functional determinant and some generalizations need to be used [25l - [27| . 
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In what follows, we will assume that the logarithmic term in the heat kernel asymptotic expansion does not appear, 
this will allow a standard definition of the functional determinant (T(j |. From a more formal point of view, as said, 
this means that C(s) is regular at s = 0, which is the case if Cyy{s) is regular at s = ~l/2- As we will show in the 
next section, this this means that [Qp]^^^ vanishes identically, since it is proportional to the residue of at the point 
s = -1/2. 

With the previous remarks in mind, from (1351) and (j36|) one can easily find that 

C(o) = [gj,],^^ + [[gpi^^-gp{o)) , (37) 

C'(0) = {[GRl^g - QRiO)) - du^ inu^l^gpiu') . (38) 

The above formulas are general. In the following sections we will consider the case of Riemann caps ^ and compute 
the terms that are necessary for the evaluation of C(0) and C'(0)- 



IV. GENERAL RESULTS 



The results of the previous section offer a way to compute the quantities C(0) and C'(0) and, in turn, the functional 
determinant. Specifically, the uniform asymptotic expansion (j25p can be used to obtain the various terms in formulas 
([57]) and ([55]) and, thus, a general expression for the determinant, as we will show in the following sections. 



A. Logarithmic contribution 



Let us begin by computing the logarithmic contribution. According to the above Lemma, we need to expand the 
function Q for large values of p = w/ik [l^. The coefficient of the logarithm will then provide the term [g]igg- This 
coefRcient is, in principle, composed of two terms [GR\iog and [Gpjiog- The first step is, then, to expand ([^5)) for large 
p. From its explicit expression one gets 



_ 1 1 °° 
InP _ (cos6'o) — (1 + 2^k) hrp^ ln(27rsin6'o) + parctan(tan6'o) + 



n=l P 



ri=l 



where T{9q) are continuous functions of the parameter 0o which are computable from the expansion of t(u) in (|23p . 
and the terms an(j^) behave like m~™ with m > 1 for large u. 

It is trivial to see that only the first term of the previous expression gives rise to log terms. The remaining sum 
over the eigenvalues is not difficult to compute according to the definition of the base zeta function, leading to the 
result 



\g\, is) = lim 

Hon \ I /I „ .n 



4 s^o 



(40) 



It is clear, from this formula, that [Qp\^^^ does not vanish if Qjy is not regular at s = ~l/2- In fact, as mentioned 
before, by denoting with Res the residue, one would obtain 



[Qp\loq 



^^R-esC-r 



ResC(O) 



(41) 



which is a relation that has been shown to hold also for the generalized cone [lOl ■ 

In order to deal with a well defined functional determinant, we will assume that ResCyK(^l/2) = 0. In this way, 
by taking the limit as s — > 0, Cyvi^ ~ 1/2) is finite and C^(s) is regular, as formula (|31l) indicates. This allows us to 
obtain the result 



[QR\log ^ 



C>'(0) + 2C^ - 



[Qp],.. = 



(42) 
(43) 
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B. Polar contribution 



The next term we need for our analysis is the polar contribution to the spectral zeta- function, namely [Gp], which is 
computed from the expansion of the expression (j25p for large /ik [l3j ■ In this way, the sum over /Ltk can be performed 
and expressed in terms of the zeta- function on the base manifold . Since the analytic structure of is known, it 
will not be difficult to extract the polar part [Gp]- It is straightforward to show that the expansion of (|25p in inverse 
powers on /ik is 



InP' 



11 °° 

■ (cos^o) ^ -In t - - ln27r^k + A*kT(u) - ^ikln Mk - A^klnu + ^ 



an{iy) 



(44) 



By performing the sum over /ik in the previous expansion, from the basic definition of the base zeta function, one 
readily obtains 



.(cos(?o) ^C^{s){\nt~\n2TT) + u(^s~^^{T{u)-\nu) + ^C'^{s) 



s 

2 



n—l 



where we have used the following relation 



lim 



= -E'^(kK''lnA*k = k:^(s) 



(45) 



(46) 



J q=q 



From the expansion (|45p we only keep the terms that will contribute to the polar part. It is not difficult to see 
that the only potential terms contributing to the polar part [Gp] come from the infinite series in (|45p . In fact, by 
exploiting the Mittag-Leffler expansion (pi]) one can soon realize that possesses poles of the first order for k < d 
at s = d/2, (d - l)/2, . . . , 1/2 with residues [H, [H H!] 



ResC^(s) 



ri 



^^j- , where n = 0, 1, 2 . . . , d — 1 . 



This remark allows us to write an explicit expression for the polar part, namely 

d 



Gp{u) = ^ a„(i')ResC^ 



(47) 



(48) 



This relation can be used in order to find tJp(O) which is needed for the computation of C(0) as formula ([57| shows. 
By taking the limit as u — > one simply obtains 



a 



(49) 



n=l 



where we have used the fact that, according to the relations in (|23p . lim„_j.o '^('") — cos^o- The terms a„((T, ^o) are 
polynomials in sin(0o/2) and tr, and have the following general form 



x„(a,^^o)=C„ + |]QV'sin2' (^^) 



(50) 



Some of the first few a„(cos0o) £^re listed in the appendix. The numerical coefficients C„ can be found by taking the 
limit as 00 ~^ of the previous expression. By analyzing, in this limit, the recurrence relation (|25p and the expansion 
(1^51) one can easily realize that 



Bn+i CRi-n) 



n(n + 1) 



(51) 



where (^p (n) is the Riemann zeta function and the last equality can be proved by taking the limit as g — > of the 
functional relation [301 



CH{-n,q) = - 



Bn+ijq) 
n+1 



(52) 



valid for n > 0, with (^H{n,q) being the Hurwitz zeta function and B„{q) being the Bernoulli polynomials. 
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C. Regular contribution 

The regular term [G]reg instead obtained by isolating the coefficient of the term in the large p expansion that 
we have obtained in (1391) . One can easily see that regular part has, then, the following form 

[gieg{s) = -lin{2TTsm9o)U{s) . (53) 
From this expression we can isolate the polar and regular parts of (s) as s — ^ 0. In more details we obtain 

[GRleg = -\ln{2iTsmeo)U(Q) , (54) 

and 

[Gpleg = ■ (55) 

The other term necessary for the evaluation of C'(0) is ^_r(0) which has the form 

gn{0)= PF|limX^d(k)Aik'MnP;^^/2(coseo)| , (56) 

where PF denotes the partie finie of the expression. It is clearly not possible to find a closed expression for the above 
summation. However, some progress can be made by expressing the Legendre function in terms of a hypergeometric 
function as follows, 



1 ... / Oo\ f 1 1 , 1 . • 2 / 



P;-^V2 (cos ^o) = ^^^^^ tan- (^y j (^-- + ^ , - + ^ - + 1; sm ^- 

By substituting the previous expression in ()56p . one obtains 

5RiO,s) = lc^(s)+\ntan(^-^)u(-s-l) - PF^ d(k)A.^'' Inr(Mk) 



(57) 



PF^d(k)A*k'^ln2Fi (^-a+i,a+i,^k + l;sin2(^|^^ . (58) 



At this point we exploit the integral representation of the function \nT{fi) [TtI. |3G| 

lnr(Aik)= Uk-^jlnMk-Aik + ^ln27r+ / dt - j + ] '-^ , (59) 



to obtain 



^ d(k)/.-2^ Inr(^k) - -X'^ [s) is) In 2^ + PF A(s) , (60) 



k 

where we have introduced the function 



In order to find the finite part of A(s) one can simply subtract from its full expression (|ST|) the polar part Apoio which 
has been evaluated in (|50l) . Since A does not depend on the parameter 6*0, we obviously need to subtract from (pTj) 
only the terms in (j50p which are independent on ^o- By doing so we find the following finite expression at s = 0. 



PFA(0) = .i.„(*W + i|;MzIl)R,e.o.(=)^ 



(62) 
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For a more explicit evaluation of (|62|) it is convenient to recast it in the following form 



PFA(0) = ^d(k) 



dt 



E 



2 t e* - 1 



5:<sL!i)ppc,(|), 

n=l 



where we have used the relation 



^ResC^ [-) = 



<.(f)-PFC.G) 



and the integral representation of the inverse powers of /ik as 

1 



Tin) 



(63) 



(64) 



(65) 



In order to deal with the integral in ([63)) we will follow the methods developed in 0, [13] where a similar object has 
been studied in the context of the d-dimensional ball. At this point, it is convenient to define the following square 
root of the heat kernel 0, HH 

i^^^(t) = ^d(/ik)e-*^-. (66) 

By recalling that the spectral ("-function is obtained from the heat kernel by inverse Mellin transform as [13, [12, l33| 

1 

we can write the expression in (j63p as follows 

P-A(o^^) ^ E^r(...)c.(i±^).ir(.)c.(0-r(.-i)c.(i^) 



(68) 



where we have introduced a regularizing parameter z such that PF A(0, z) — >■ PF A(0) as z — 0, and we have introduced 
a new function defined as 0, l3lj | 



r(z) 



(69) 







The limit of the above expression as z — >■ can now be taken leading to the result 

d 



PFA(O) = 



n=l 

7 



2*(n)ResC^(|) + PF (|)] + JcV (0) + ^CV (-^) +U (-^) 



(0) - lU 









+ lim — 







2E 



— ^ResC^ 



C^(0) 



+ Cyy ( --) +r(z)CAA+i(2) 



E 



(70) 



where ^'(n) is defined as the logarithmic derivative of the gamma function. At this point, it is convenient to utilize 
the series 



1 1 fi 



5* - 1 t 2 ^71 

n— 1 



1/2 

and the asymptotic expansion of the heat kernel K (t) for small values of t as 



(71) 



(72) 
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where 



4 = 2r(d- fc)ResC^ 



d-k 



d-k 



for yfc = 0, 1, • • • , d - 1, yfc = d - (2Z + 1) with leN+, and 
for fc > rf and k € With the help of the above results, it is not very difficult to prove that 0] 



-C^(0)-2^ 



Res 



(73) 



(74) 



(75) 



The last expression ensures that the limit as z — > in ()70p vanishes identically, as expected, yielding the final result 



PFA(O) = c;,+i(0) + -CV(0) + ^CV +Ca^ +2E 



E^I^-C^g), (76) 
,fe=i J 



where in the formula ([75]) we have used the fact that [3C 



^'(?i) 



n — 1 ^ 



fc=i 



(77) 



In order to deal, now, with the last term in expression (|58l) we use an expedite method relying on the Abel-Plana 
summation formula fl§L [l3 | . Let us define the function 



i^(x) = d{x.) \n2F1 ( -o- + i (7 + + l;sin^ 



(78) 



From the series representation of the hypergeometric function and from the fact that (i(x) is polynomial in x, one can 
see that our function 7(x) uniformly satisfies the condition limj,^o e^^'^'^' \"f{x + iy)\ = 0. This allows us to apply the 
Abel-Plana formula as follows 13, 17] 



PF ^ d(k)Aik'^ In 2F1 (^-a + i a + i, /ik + 1; sin^ (^|^ ^ 



PF 



m: 



2s 



dx + i 



(p(ix) + w(—ix) , 1 

—dx + -ip{0), 



(79) 



where he have set s = in the last two terms because they are finite for s — >■ 0. The finite part of the first integral 
in (|79l) can be obtained in the same way that was described for the function ^{s) in (pT|) . By subtracting from the 
s-dependent integral in (j79p only the polar contributions in (j50p depending on 9q one obtains 



PF 



//2s 
Mx 



(80) 



At this point, the above integral can be analyzed further, by utilizing repeated integration by parts, once the base 
manifold .yV is specified and the degeneracy of the hyperspherical harmonics on ^/K is known (see e.g. [13, 113]) • 
By collecting the relevant results obtained above, we can finally write an explicit expression for Gr{0) as follows 



PF 



In tan 

(p{x 



ic^(0)ln2.-CVi(0)-2E^^ ResC^^Q) 

n=l Lk=l J 



Mi 



2s 



dx + i 



(p(ix) + ip(—ix) , 1 



- 1 



(81) 



The last term that is necessary in order to evaluate C'(0) is the integral in (p8)) . which represents the non-local 
contribution to C'(0)- By recalling the polar terms already found in P5|) . we can write that 



du^lnu^— ^gp(ii2) = VRcsC^ (-) / du'^lw, 

^ 2 / Jq 

The functions a„(i^(w)) can be computed, as many as needed, from the cumulant expansion (|26p . 



(82) 
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D. C(0) and C'(0) 

We are now able to write down an expression for C(0) and C'(0) by using the results obtained in the previous 
sections. By recalling the expression (|57)) and the results (H^ . and (1551) "^"^ have the following formula for the 
value of the zeta function on 17 at s = 



C(o) = - \u 



a 

^ a„(cr, 6'o)ResC^ 



(83) 



where the polynomials an{cr,Oo) are given in ()50p . In order to write an expression for the derivative of C(s) at s = 0, 
we need the formula psp and the results obtained in Sec. IIV CI The final expression is slightly more cumbersome 
and reads 



C'(o) 



ln(sin6'o) + 



1 



ln(l+cos^?o) + 0+i(0) 



2E 



Cni-n) 



n=l 



2^^ ' .L e2-- 



E ^ Res (^) - ^ Res (f ) / In 

fc=l J ri=l 

da; - hm <^ / -^dx > ResC^ / . 



^'^sin^'f ^ 



(84) 



where the function Lp{x) has been defined in (1781) and the expression in the last line is finite, by construction, as s ^ 0. 

Although the above expressions are quite implicit, they are very general and valid for arbitrary dimension and for 
any smooth and compact base manifold . We would like to stress, here, that more explicit formulas can be obtained 
once the base manifold is specified. In the next section we will consider a particular case in which the base manifold 
is a d-dimensional ball. This case is of special interest because the spectral zeta function on d-dimensional ball can 
be explicitly evaluated in terms of the Barnes zeta function [l3| • 

V. A PARTICULAR CASE: d-DIMENSIONAL BALL AS BASE MANIFOLD 



In this section we will consider the case in which the base manifold is a d-dimensional ball. In this particular 
situation the eigenvalues of the Laplacian on ^ are known to be 



d-l 



with fc > 0, and the eigenfunctions are hyperspherical harmonics with degeneracy 

{k + d-2)\ 



d(k) = (2fc + (i- 1)- 



k\{d-l) 



(85) 



(86) 



The explicit knowledge of the eigenvalues /^k and their degeneracy allows us to write, by using the definition (jl8p . the 
zeta function on ,yV as follows 



. , ^ Y^/o, , + f, d-l 

Q^{s) = 2_^{2k + d- D-^TT^ rri- I k 



fc=0 



kl{d-l)\ 



-2s 



(87) 



It is straightforward to show, with the use of algebraic manipulations on the factorials, that the above zeta function 
can be written in terms of a sum of Barnes zeta functions la, [lOl 



. , ^ d+l\ ^ d-l 
U{s) = Cb 2s, —- + Cb 2s, —- 



where CBis,a) is defined as |35l. ISC 



1 

(a + m • r)^ 



(89) 
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valid for 5ft(s) > d where m and r are rf-dimensional vectors, and where the notatfon o.\1) = Cb{s, a) has been used. 
In order to extract the information we need from the Barnes zeta function we will utilize its integral representation 
as follows O,!!! 

Cs(s,a) = — dy , (90) 

27r 2'*smh''(|) 

where L represents the Hankel contour. With the help of the above representation and the relation (|88l) we can write 
an expression for C^i^) ^ 

^ ' 2^ Jl ^ 2'i-isinh'*(f) ^ ' 

At this point it is suitable to make a change of variables, y/2 — > y, to recast the previous formula in the form [lol. Isij 



^ I c\ I 1 ^2s-2 ^r(2 — 25) 2s+l-rf \ " f 



dyy^^-^-^+- ^ (92) 



where the coefficients dI'^ '^'^ can be easily computed, by equating like powers of y, from the formula 



sinhy 



d—l 00 



= J2Dl^-^^y-. (93) 



The integral representation (|92p is particularly useful for computing the residue of at different values of s. 
Specifically, we need the residue of C^i^) the points s = m/2 with m being a positive integer. It is not very 
difficult to see that the integral in (|9 2l) vanishes unless v ~ d — m where the integrand has a simple pole. By simply 
using the residue theorem we get [31j 



Res( 



/ TTl\ '2^m — dj^i^ 1) 

'^-^(y) " (d-l)(TO-2)!(d-m)! ' ^^^^ 

valid for m > 2 and d > m. Moreover, we will need the values of at the points s = and s = —1/2. From the 
expression one easily gets 

ryl-djjid-l) 

(-5) - (i -mlxf ^ 

These last formulas immediately give an expression for C(0) for arbitrary dimension d. More precisely, by recalling 
(l50l) and (|83)) . one obtains 



and 



C(o) 



+ 1 



2rf+i(d- l)d! 
1 



1 T^U{-n)Dtn^ 
2\d - 1) ^ n(n - 2)!(d - n)! 



1 r^n T^{<i--'^) n 



(97) 



Let us, now, turn our attention to the term C^/k^i(s) that appears in the general result (j84p . In the case of a 
d-dimensional ball as base manifold we have that |10i] 

C^+i(.)=^e(fc) fc + ^ , (98) 
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where we have introduced the coefRcients 



3(fc) = {2k + d) 



{k + d-iy. 
kidi ■ 



By writing 0, |3l| 



e(fc) = ^e„(d)(fc 



Q = 



d+1 



(99) 



which defines ea{d), we have that 



Q = 



<i+ 1 



(100) 



where the coefScients (d) depend on the dimension d and are determined from the equation (j99]) , and Ch represents 
the Hurwitz ^-function. In particular, for its derivative at s = we obtain 7„ ,31.] 



(101) 



This expression, involving the derivatives of the Hurwitz zeta function, can be rewritten in terms of Riemann zeta 
function and its first derivative. In order to do so, we must distinguish between even and odd dimensional cases. For 
even dimensions, namely when d — 2q with q G N^, we get (3l| 



d+1 



= CH[-a,q+-] = -7^ 



In 2 



n=l 

2q~l q-1 

+ — In 71- In 

n— 1 n—1 

while for odd dimensions, i.e. when d — 2q-\-\ with q G N^, we obtain 



+ (2-- - l)Ck(-") 



(102) 



d+1 



9-1 



(103) 



n=l 



We would like to mention that the derivatives of the Riemann zeta function appearing in (|102p and (jl03l) can be 
further simplified by differentiating the reflection formula 



Cfl(l - 5) = 2(2^)^™ cos (^) Cfl(^*) . 

In fact, for a = 2p where p G N+, it is not very difficult to get js^, 

Cr{-2p) = (-l)f ^(2^)-2f-ir(2p + l)C;^(2p + 1) , 
and for a = 2p + 1 where p E N+ , one obtains 

H- 1) ^ f [.(2.) - m 2.] + ^^«^C^(2.) . 



(104) 



(105) 



(106) 



The next term that can be more explicitly evaluated is the integral in ([5D|) . By recalling the expression for the 
degeneracy of the eigenvalues (|86)) we can rewrite the function (|78p as follows 



,(a;) =x'^^e„(d)a;"-'^ 1 



d-1 
2x 



In 2F1 ( -o- + -, (7 + -, /ix + 1; sm 



(107) 
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This formula allows us to write 



ip{x) 

,2s 



dx 



where we have performed a change of variables to 1/a; and we have defined 



0(x) = ^e„(d)a;'^-" 1 



d-1 



X \n2F1 -cr + -,CT 



1 11 



2x 



1; sin 



■ 2 n^o 



(108) 



(109) 



which is an analytic function at x = 0. It is convenient, at this point, to split the integral (|108p into three pieces 
[11 [3 as 



7(2;) 



Mi 



2s 



dx 



dx X 



x) + dxx~'^~'^(j){x) 



dxx~'^~^(j){x) 



(110) 



where we have set s = in the last two integrals because they yield a finite result in that limit. The first two terms 
can be integrated by parts d + 2 times to get [H [l3 | 



7(2;) 



Mi 



2s 



dx =^ 



2s - d - 1 

+ •••-(-1)'^ 

+ (-1)' 



{x)- 



d(j){x) 



(2s-d- l)(2s-d) dx 



„2s 



d'^+^^{x) 



{2s-d-l){2s~d)---2s dxd+^ 



dxx^^''^'^^^^ 



(2s -d- l)(2s-fi)---2s 



dx<i+2 



'd-i 



(-d-1) 
•-(-1)' 



0(x) 



d(j){x) 



{-d-l){~d) dx 

\nx d'^+^(f>{x) 
{-d - l){-d) ■ ■ ■ (-1) dx<^+^ 



J 1 



(-1)"^ 



(-d-l)(-d)...(-l) A 



dx In 



^rf^+20(x) 

' dxd+-^ 



dx x 



-d-2 , 



(111) 



We can now perform the analytic continuation to s = of the expression above to get the compact result 



m: 



2s 



dx = 



1 d'^+^(j){x) 



2s(d+l)! dx-i+^ 
1 



{d+l)\ 



dx In X 



x=0 

d'^+'^(j){x 
dx'i+^ 



id 



dx'^+^ 



x=0 



dxx-'^-^(f>{x) 



(112) 



where the dependence on the intermediate cutoff a; = 1 has disappeared as expected. From this last formula, by 
recalling (jSO]) . one can easily find the relations 



and 



PF 



M: 



n=l 



^dx = - 

2- [d 



1=1 

fc^O 



1 d'^+^(t>{x) 



2(d+l)! dx-i+i 



(113) 



x=0 



+ 1)1 [tod-k + l] 



dxx-'^-^cjjix) . 



d'^+^(f){x) 



x=0 



(d +1)1 Jo 



dx In : 



d'^+^(f>{x) 
' dx<^+^ 



(114) 
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A. Specific dimensions 



In this section we will list the results for ^(0) and C'(0) for a dimensional ball as the base manifold of O in a 
few lower-dimensional cases. Let us start with d = 2, namely a 2-dimensional base, meaning that is of dimension 
D = 3. In this case we have 



and also 



Ci.=3(0) = - 3^ (-1 + sin' ^0 , 



J du^ In -^-^a2{i^{u^)) — i 



(fi{ix) + ip{—ix) 



- 1 



dx + 



x=0 



dxln. 



d'^(j){x) 



POO 

+ / dxx~'^(j){x) . 



For d = 3, i.e. D = A, we obtain 

1 



1 



180 8 



(-l + 4c7^)sin2 



64 



(-l + 4c7^) (-25 + 4c7^)si 



sm 



- ^(-l + 4a^)(-25 + 4a^)sin«(|) 



and 



1 



Cd=4(0) = ^ - ^ ln(sin 9o) + ^ ln(l + cos Oo) + Jck(-l) " ICni-^) + ^C^? (-3) - ^^(0) 



25 




288 


2 






1 


d^6(x 


— / dxinx V- 


24 Jo 


dx° 



du^ Inu^— raafi/fu^)) — i 
dw^ 



(p(ix) + (p{—ix) 



- 1 



dx 



+ 777/ c^a;lnx — h 

24 io dx^ 

For d = 4, or equivalcntly D — 5, one gets 
17 



dxx ^(j){x) 



+ ^(-l + 4»')(-13 + ta2)smi 



and, finally, 

Cd=5(0) 



5^ ^o) + ^ ln(l + cos .„) - y,i-,) + lc;.(-3) - ^C;.(-2) - ^^Cni-D 

1 _ 137 d'^<t){x) 
2*^^ ^ ~ 7200 dx^ 



'f 

Jo 



if{ix) + if{—ix) 

g27rx _ 1 



,T=0 



dx 



1 Z""" d 

24 y '^"^^'^"^rf^? [«2(!^(w')) - 4a4(!^(w'))] 



(115) 



(116) 



(117) 



(118) 



(119) 



(120) 



It is quite straightforward to see that the above results reproduce the known results for the ball once the limit for 
small 9q is taken. Moreover, for conformally invariant scalar fields it is not difficult to realize that the parameter a 
takes the value 1/2. As one can easily see from the above formulas, when this case is considered our results for (^(0) 
reduces to the ones obtained in flat space as it should be expected. We would like to stress, here, that although the 
integrals appearing in the above formulas for C'(0) are not explicitly computable in the general case, one could start 
immediately a numerical analysis obtaining, in this way, very specific results. 
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VI. CONCLUDING REMARKS 



In this paper we have presented a technique for the evaluation of the zeta determinant for the massive Laplacian 
on a Riemann cap (or spherical suspension). Riemann caps are generalization of deformed spheres in the sense that 
the base manifold is compact and boundaryless, but otherwise general, and present a singular Riemannian structure 



analogous to that of generalized cones. In the spirit of ref. [18[, we expressed the functional determinant in terms of 
the base zeta function and found some general formulas for the analytical continuation of the zeta function and its 
derivative at s = 0, which may be used to obtain the determinant and, in a more physical context, the effective action 
for a quantum field propagating on the above manifold. Finally, we specialized our general formulas to the case of 
a ball as a base manifold and presented the, still general, results in such case. Explicit formulas are finally given in 
D = 3, 4, and 5 dimensions recovering previously obtained results. 

Aside of presenting results which are more general than those available in the literature, one of the goals of our 
work was to illustrate a method to compute the zeta determinant alternative to that of ref. [id ] , which is the one now 
commonly used in various situations. The basic idea behind the method of ref. [lo| is to express the zeta function 
as an integral in a way similar to what we did here, and then add and subtract the asymptotic expansion of the 
integrand. This method of asymptotic subtraction allows a direct way of isolating the divergences and thus regularize 
the functional determinant by direct computation of the integral of the asymptotic part. In general, this last step 
depends on the actual behavior of the egeinfunctions and it is non-trivial. In fact, in the case of Riemann caps analyzed 
here, it does not seem to be the most convenient way to proceed. Alternatively, the method we have used here is based 
on an adaptation of a Lemma of analytic continuation, expressed through formulas ([33l) - p4|) . that relates the zeta 
function and its derivative at s = to the coefficients of the (uniform) asymptotic expansion of the eigenfunctions. 
It goes without saying that it can also be used in the well known cases of spheres and balls, reproducing previous 
results. 

One situation of physical interest, where the present technique may be used, is in the context of higher codimension 
brane- world models. In these models cone- type singularities arise and one standard way of regularization is to cut the 
tip of the cone, thus removing the singular behavior [SOi - fi^ . In fact, the case of conic manifolds has received some 
recent attention, and, for instance, ref. [23| describes some unusual and pathological behavior of the zeta function that 
cannot be easily recovered by simply taking the limit of cone regularized in the above way. One possible alternative 
is to smoothen the cone by rounding it with a spherical cap. In this way, recuperating the conical structure is rather 
transparent and can easily done by taking the limit of small angles, where the metric of the Riemann cap reduces to 
that of the cone. 

Another advantage of the method used here is that the functional determinant is expressed in terms of the coefficients 
of the asymptotic expansion ([33]) making the method suitable for numerical implementations. This may be particularly 
useful in computation of the Casimir energy. In fact, the method utilized here can be easily generalized to consider 
the analytic continuation of the zeta function to other points in the complex plane allowing to: obtain the analytic 
continuation of the zeta to other values of s and thus constrain further heat-kernel coefficients; develop a different 
way to compute the Casimir energy, which requires the value of the zeta determinant at s = —1/2. 

The above issues are currently being considered and we hope to report on our progress soon. 
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Appendix A: Polynomials a„(cos6o) 

For the reader's convenience, we list, here, the polynomials a„(cos0o) up to the fourth order. By utilizing the 
recurrence relations (|25l) and (1^51) and by subsequently taking the limit as u — one obtains 



1 1 



«i(^,^o) = -^ + ^(l-4a^)sin^(^) , (Al) 



a,{a,9o) = -i(l-4a^)sin2(|)+i(l-4a^)sin^(^ 



(A2) 
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03(0-, ^o) 



^ + i(l-4a^)sin^ 



2 / t'o 



4 / ''o 



— (25 - 1040-2 + 16(7*) sin , „ 
32 ^ ^ V 2 



+ 1(25-I04a2 + I6a4)sin<^(|) , 



(A3) 



a^a, ^o) = (1 - 4(7^) sin^ (^|) + ^ (^^ - 640^ + 16a*) sin^ 

- ^ (13 - 56a2 + 16a*) sin^ f + 1 (13 - 56a2 + 16a*) sin^ f ^ 
4 \^ 2 y 8 \ 2 



(A4) 
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